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ABSTRACT 

The nature of the broad-band spectra of supermassive accreting black holes in active galactic 
nuclei (AGNs) is still unknown. The hard X-ray spectra of Seyferts as well as of Galactic 
stellar-mass black holes (GBHs) are well represented by thermal Comptonization, but the 
origin of the seed photons is less certain. The MeV tails observed in GBHs provide evidence 
in favour of non-thermal electron tails and it is possible that such electrons are also present in 
the X-ray emitting regions of AGNs. 

Using simulations with the kinetic code that self-consistently models electron and pho- 
ton distributions, we find that the power-law-like X-ray spectra in AGNs can be explained 
in terms of the synchrotron self-Compton radiation of hybrid thermal/non-thermal electrons, 
similarly to the hard/low state of GBHs. Under a very broad range of parameters the model 
predicts a rather narrow distribution of photon spectral slopes consistent with that observed 
from LlNERs and Seyferts at luminosities less than 3 per cent of the Eddington luminosity. 
The entire infrared to X-ray spectrum of these objects can be described in terms of our model, 
suggesting a tight correlation between the two energy bands. We show that the recently found 
correlation between slope and the Eddington ratio at higher luminosities can be described 
by the increasing fraction of disc photons in the emitting region, which may be associated 
with the decreasing inner radius of the optically thick accretion disc. The increasing flux of 
soft photons is also responsible for the transformation of the electron distribution from nearly 
thermal to almost completely non-thermal. The softer X-ray spectra observed in narrow-line 
Seyfert galaxies may correspond to non-thermal Comptonization of the disc photons, predict- 
ing that no cutoff should be observed up to MeV energies in these sources, similarly to the 
soft-state GBHs. 

Key words: radiation mechanisms: non-thermal - accretion, accretion discs - galaxies: 
Seyfert - X-rays: galaxies 



1 INTRODUCTION 

Accreting radio-quiet supermassive black holes (SMBHs) resid- 
ing in the centres of quasars, Seyfert galaxies, narrow-line Seyfert 
1 galaxies (NLSyl) and some fraction of low-ionization nuclear 
emission-line regions (LINERs) in many respects are analogous 
to the Galactic stellar-mass black holes (GBHs) in X-ray binaries. 
In the X-ray/soft 7-ray band, the spectra of Seyferts can be repre- 
sented by a sum of a power-law-like continuum which cuts off at 
a few hundred keV and a reflection component with the iron flu- 
orescent Ka line believed to be produced by reprocessing of the 
intrinsic power-law by cold opaque matter, probably the accretion 
disc (Nandra & Pounds 1994). The 2-10 keV intrinsic spectral en- 
ergy slope (defined as Fe oc E'") of the power-law a ~ 0.9- 
1.0, which is ubiquitously found in Seyferts, is somewhat larger 
than what is measured in GBHs in the ir hard state a ~ 0.6-0.8 
(e.g. lZdziarski. Lubinski & SmitrilT999l) . NLSyl having softer X- 



ray spectra than Seyferts, probably represent a state with higher ac- 
cretion ra te (in Eddington units) similar to the soft/very high state 
of GBHs jPounds. Done & Osbornj|l995h . 

The seemingly similar X-ray slopes of Seyferts triggered 
the efforts to find a physical model which would explain such 
spectral stabilit y. The non-thermal models became popular (see 
^Sven ssonlll994 for a review). These invoke injection of high- 
energy leptons into the system with subsequent Compton cool- 
ing by accretion disc photons and photon-photon pair produc- 
tions initiating pair casca des. Saturated cascades produce intrin- 
sic sp ectra with a ~ 0.9 jZdziarski & Lightmanlll98^ : ISvenssonl 
Il987l) . which becomes consistent with the X-ray spectra of Seyferts 
after accounting for th e hardening due to Compton reflection 
( lGeorge&Fabianlll99ll) . The pair-cascade models, on the other 
hand, predict a strong tail above 300 keV and an annihilation 
line which have never been observed. The CGROIOSSE ob- 
servations of the brightest Seyfert 1 galaxy NGC 4151 con- 
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strain the fraction of energy going to non-thermal injection to 
be less than 50%, while the rest of the power going to ther- 
mal heating (Zdziarski, Liehtman & Mac iolek- Niedzwieckill 19931 : 
Izdziarski. Jo hnson & Magdziarz 1996). The existing upper limits 
on the average flux of Seyferts above 100 keV a re compatible 
with the presence of wea k non-thermal tails (e.g. iGondek et al.l 
Johnson et al. 1997h . The more recent INTEGRAL observa- 



Il99d 



tions jLubiiiski et al, 



2010 ") have not improved those constraints. 



Non-detectable high-energy tail in Seyfert galaxies can also 
be interpreted in terms of pure thermal models, where the power 
is equally shared among all the thermal particles. The spectra in- 
deed can be well described by Comptoni zation on thermal 50- 
150 keV electro ns (see e.g. Z dziarski etaPl 19971 : |Poutaneniri998l : 
IZdziarski|[T999h . The stability of the spectral slopes can be inter- 
preted as the evidence of the radiative feedback between the hot 
X-ray emitting plasma and the cool accretion disc, which is as- 
sumed to be a sole source of the seed photons for Comptoniza- 
tion. The slab-corona mod el, where the hot plasma san dwiches 
the accretion disc (see e.g. iHaardt & Maraschi|[l99ll Il993i) . pre- 
dicts too soft spectra due t o large flux of the reprocessed UV 
photons i lStem et alJ [l99l . IaM zac. Dumont & MoucheJ j2005h 
showed that the hard spectrum can be achieved if one assumes large 
ionization parameter of the disc, but in this case the model fails on 
the predicted reflection properties. Localized a ctive r egions atop 
of a cool disc llJaa rdt. Maraschi & Ghisellinii jl994 IStem et al] 
1 19951 : 1 Poutanen & Svenssonlll9 96: Svensson 1996) are more pho- 
ton starved and can produce much harder spectra in agreement with 
observations. However, in this case it is difficult to understand why 
there is a preferential slope, as it is a strong function of the separa- 
tion of the active region from the disc. 

The detection of the MeV tails in the spectra of the GBH 
Cyg X-1 in both hard and soft states l lMcConnell et alj |2003) 
imply the presence of a significant non-thermal component in 
the electron distribution. The power-law looking spectra extend- 
ing to the MeV ene rgies in other soft-state GBHs (Grove et al.l 
ri998i : rzdziarski et al.ll200]]) ar e consistent with being produced by 
non-t hermal Comptonization l lPoutanen & Coppl 1 19981 : IPoutaneiil 
1 19981) . By analogy, the electrons in Seyferts also might have a 
non-thermal population. More probably in both types of sources, 
the electron distribution is hybrid in bo t h hard and soft states 
jPoutanen & Coppilll998l : |Poutanen|[l993 : ICoppill 19991) , with the 
non-thermal fraction increasing for softer spectra. However, non- 
thermal particles could be spatially separated from the thermal 
ones. 

The presence of the non-thermal particles, even if they are 
not energetically dominant, has a strong impact on the emitted 
spectrum, because the synchrotron emission can increase by or- 
ders of magnitude ev en if only 1 per cent of electr on energy is 
in the power-law tail dWardzinski & Zdziarskill200lh . As most of 
the synchrotron emission is self-absorbed, this process plays an 
important role in shaping the electron distribution by thermaliz- 
ing particles via emission and absorption of synchrotron photons, 
the s o called synchrotron boiler (Ghisellini, Guilbert & Svensson 
[T98g : iGhisellini. Haardt & Svenssorjl 19981) . 

It was recently shown that under the conditions of GBHs, 
Coulomb collisions and the synchrotron boiler efficiently thermal- 
ize electrons producing thermal population at low energies even if 
the electrons were originally non-thermal (Vu rm & Poutanen 2008; 
IPoutanen & VumJ 20091 : [Malzac & Belmonfcoogi) . The hard spec- 
trum of GBHs can be fully accounted for by the synchrotron self- 
Compton (SSC) mechanism in the resulting hybrid electrons. In 
this picture the disc photons are not required to interact with the hot 



X-ray emitting plasma at all. This mechanism also gives a rather 
stable X-ray slopes for a large range of parameters, relieving the 
need for the feedback between the disc and the active region. 

As many of the radiative processes depend only on the com- 
pactness of the source, but not on its luminosity or size separately, 
one may think that the spectra of Seyferts can also be described 
in terms of hybrid Comptonization models. However, not all pro- 
cesses can be scaled away. The aim of this paper is to study in 
details the spectral formation in relativistic hybrid plasmas in the 
vicinity of a SMBH. One of our goals is to learn how the results 
scale with the mass of the compact object. We also test the role of 
bremsstrahlung in SSC models and, finally, we compare the model 
with the data on Seyferts and NLSyls. 



2 MODEL DESCRIPTION 
2.1 Geometry 

It is commonly accepted that the main process responsible for 
X-ray spectra of radio-quiet active galactic nuclei (AGNs) is 
(nearly) the thermal Comptonization of soft seed photons by en- 
ergetic electrons. The nature of the seed photons is less under- 
stoo d. The two main candidat es are the standard cool accretion 
disc dShakura & Sunvaevll973h and the synchrotron emission from 
the hot electron populations. We model both types of seed photon 
sources. 

The geometry of the X-ray emitting region also remains un- 
clear. It can be the inner hot part of the accretion flow which 
ca n be radiatively inefficient as well a s efficient (s ee reviews 



corona models are also popular ( Galeev, Rosner & Vaiana 19791: 


Haardt & Maraschi 1991. 


1993; Haardt etal. 


1994; IStem et all 


I995I: IPoutanen & Svensson 1996: Poutanen 


199^. The slab 


(plane-parallel) corona jHaardt & Maraschill99] 


,1 1993b is unlikelv 



to be realized, as it produces Comptonization spectra which are 
too soft compared to the observed a ~ 1 (Stern et al. 1995). The 
corona can also be patchy and consist of iso lated active r egions 
atop a cool accre tion disc ( lGaleevetal.lll979l ; lllaardt et alii 1994 
IStern et alJ[T99l) . This kind of geometry can produce spectra of 
almost arbitrary hardness if only soft photon reprocessed in the 
disc are considered. The corona does not need to be static, but 
ca n have substantial bu lk velocity away from the disc or towards 
it dBeloborodovll 19993 ). In this work we consider spherically sym- 
metric region, which can be either the inner hot part of accretion 
flow or a blob of gas above the accretion disc. No anisotropy in 
geometry or seed photons is taken into account. 



2.2 Main parameters 

In most of the simulations we assumed that the active region has a 
size of R— 10-Rs, where Rs — 2GMbh/c^ is the Schwarzschild 
radius of the black hole of mass A/bh, as one can expect most 
of the energy to be liberated within a region of a corresponding 
size. It is convenient to introduce a compactness parameter, which 
is independent of the black hole mass: 
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where LEdd = 47r^eGMBHmpc/crT = 1.3 x IO^'^Mbh/M© erg 
s~^, (tt is the Thomson cross-section, fic = 2/(1 + X), and X is 
the hydrogen fraction. 
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The importance of synchrotron processes is determined by 



C/b-R c/Linj, where (7b 



magnetization parameter 77B 

BV(87r) is the magnetic energy density and -Linj ~ "3"^ cf^rad 
(Kad is the radiation energy density). The equipartition of mag- 
netic and radiation energy densities occurs at t/b ~ 0.25. The mag- 
netic field can be expressed in terms of adopted parameters as 
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where M7 = Mbh/IO^M©. 

We assume that the released energy is given to the elec- 
trons. The transfer of energy is modelled as a power-law injection: 
dNc/{dt dj) oc 7~'"'"J extending from low-energy Lorentz fac- 
tor 7ic = 1.0 to high-energy cut-off at 7ho ~ W^. The Thomson 
optical depth of the electrons is chosen at the beginning and fixed 
during each simulation, i.e. the number of particles is conserved. 
The injected electrons increase the optical depth; therefore, to keep 
it constant throughout the simulations, the same amount of elec- 
trons is removed from the system. This action simulates the effect 
of re-acceleration. The total luminosity of the system L is equal to 
the net power Linj = ^-R^^o(7inj — 7cq)'"cC^, where j^^^- and 
7gq are the mean Lorentz factors of the injected and the equilibrium 
distributions. 

We also consider an additional source of soft photons from 
the accretion disc in the form of the blackbody radiation injected 
homogeneously to the system. The parameters that describe them 
are the temperature Td and the ratio of the disc luminosity to the 
injected power / — I/disc/^inj- 

Thus the main parameters of the model are the mass of the 
compact object, A/bh, size of the active region, R/Rs, magneti- 
zation parameter, rjB, Eddington ratio, L/L^dd, ratio / of the disc 
seed photon luminosity to the injected luminosity, injection slope 
Finj and the Thomson optical depth r = Ra-rric. 



2.3 Characteristic time-scales 

There are three time-scales important in our simulations: light- 
crossing time tic, cooling time tcooi and re-acceleration time tm- 

The re-acceleration time is calculated for a given optical depth 
and required injection energy as 



^ ^ 47r R r7i„j 

3 C /jni 



(3) 



corresponding to the average time of the removal of the electrons 
from the system (as modelled in the code). For large optical depth, 
the re-acceleration time is larger than the light-crossing time and 
thus it can be associated with the escape time. In the case of low 
optical depth, tra < Uc, and it cannot be thought of as the escape 
time itself, but interpreted as the time between re-accelerations. 
The synchrotron cooling time can be calculated as 
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where 7s is the synchrotron cooling rate. Similarly, one 
can determine the coolin g time by Compton scattering (e.g. 
iRvbicki & LightmadI 19791) 
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The ratio of the two time-scales is roughly 
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For most of the simulations, the synchrotron cooling time is com- 
parable to (or more than) the Compton cooling time. 

The e-e Coulomb losses a t the equilibrium can be estimated as 
(e.g. lNavakshin & Meliij|l998h 
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where In A is the Coulomb logarithm (we take InA = 16). The 
bremsstrahlung cooling time is larger (approximately by the inverse 
of the fine structure constant l/ofs ~ 137). At higher energies, the 
cooling is determined by radiative processes, while for lower ener- 
gies non-radiative Coulomb collisions are dominant. The relevant 
scaling is 
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The ratio of the synchrotron cooling and re-acceleration times 
can be approximated as (7i^j ~ 2 in our simulations): 

^cool.s ^ ^'"^ /Q\ 

tra rriB 7 + 1 

and the ratio of the Coulomb cooling to re-acceleration time can be 
written as 
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As long as the optical depth is not too low, the average cooling time 
is less than the re-acceleration time. This condition is satisfied for 
the majority of simulations presented in this paper. 

2.4 Radiative processes and numerical simulations 

The gravitational energy released as matter goes down the poten- 
tial well of the black hole heats the protons and can be transferred 
to electrons via various mechanisms, for instance. Coulomb col- 
lisions with protons, magnetic reconnection, collective plasma ef- 
fects or shock acceleration. The accelerated electrons subsequently 
cool, producing a radiation field in the region. The presence of the 
latter strongly affects the electron distribution. Due to acceleration 
processes operating in the system, the electron distribution might 
not be Maxwellian. Even a small non-thermal tail (with the energy 
content of only 1 per cent) can significantly increase the net (after 
accou nting for self-absorption) sync hrotron emission up to ~ 10^ 
times jWardziiiski & Zdziarskil200ll) and therefore affect the spec- 
trum of the Comptonized component. 

For a given particle distribution (Maxwellian, power-law or 
hybrid) one can easily calculate the emission spectrum, but the for- 
mer in turn depends on the radiation field. Thus, to find an equilib- 
rium spectrum, which is self-consistent with the electron distribu- 
tion, one has to solve coupled time-dependent kinetic equations for 
all the particle species present in the system. 

Formal equilibrium solution of the kinetic equation for elec- 
trons, when cooling domina tes over the escape, can be written as 
( iBlumenthal & Gouldll 19701) 



rfiVcq 
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(11) 



where dNaq/d"f represents the equilibrium particle distribution and 
7 is the sum of all possible cooling rates. In the radiative (Compton 
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and synchrotron) cooling dominated regime, 7 cx 7" 
electrons exchange energies predominantly by Coulomb collisions, 
then 7 oc 7" . Substituting the cooling rates into equation dl lb and 
integrating (hereinafter we consider Finj ^ 2), one gets the equi- 
librium electron distribution in the form dNcq/d-y oc 7"" with 



Tinj + 1, 7 oc 7^ (Compton, synchrotron) 
Tinj - 1, 70c 7° (Coulomb). 



(12) 



In a general case, it is necessary to account for particle heating 
and thermalization; therefore, it is impossible to find an analytical 
solution; the problem has to be solved numerically. 

We consider electrons, positrons and photons. Neutrality of 
the plasma is provided by protons, which are assumed to be at 
rest and participate only in electron-proton bremsstrahlung. The 
following elementary processes are taken into account. The par- 
ticles interact with each other via Coulomb (M0ller/Bha-bha) col- 
lisions. The cooling/heating of the electrons (and positrons) occurs 
by emitting and absorbing synchrotron radiation, bremsstrahlung 
emission (electron-proton, electron-electron, positron-positron, 
electron-positron) and Compton scattering. The seed soft photons 
for Compton scattering can be provided by the synchrotron as well 
as by the accretion disc. Photon-photon pair production and pair 
annihilation are also taken into account; however, these processes 
play a minor role for considered parameters. 

We assume particle and photon distributions to be homoge- 
neous and isotropic. The calculations are done in a one-zone geom- 
etry with tangled magnetic field. The photon escape from the re- 
gion is modelled by the escape probability formalism. The escaping 
photon distribution can be described by the differential luminosity 
Le = dL/dE or differential compactness Ie = {h-nj/L) Le- 

To obtain the equilibrium spectrum, a system of time- 
dependent kinetic equations for the distribution functions of 
photons, electrons a nd positrons is solved us ing the code de- 
scribed in details by IVurm & PoutaneiJ ( |2009|) . extended to in- 
clude bremsstrahlung pr ocesses. This code is similar to the previ- 
ous codes of tha t sort by Coppi ' (1992, 1999*), but includes the syn- 
chrotron boiler jOhisellini et al.. 1988. 199 81) and treats Coulomb 
thermalization more accurately due to a grid in electron momenta 
extended to low values. A code very similar to ours was recently 
developed bv iBelmont. Malzac & Marcowith C2008.) . 



3 RESULTS 

In this section we present the detailed investigation of the features 
of equilibrium photon and electron distributions arising from varia- 
tions of different parameters. We consider pure SSC models in Sec- 
tion |3T| In Section [T2I we investigate the role of bremsstrahlung 
in SSC models. Spectral formation with contribution of the accre- 
tion disc photons is described in Section [33] Specific parameters 
of each simulation can be found in Table [T] 



3.1 SSC models 

Let us first consider pure SSC models putting f — 0, i.e. with no 
contribution from the accretion disc. We study how variations of the 
parameters influence the spectral energy distribution (particularly, 
the slope a in the range 2-10 keV) and the electron distribution. 
For the fiducial parameters, we take Mbh ~ IO^Mq, Finj — 3.0, 
r = 1.0, 7)B = 0.1, L/I/Edd = 10"^ (see Fig.^ sofid fine). The 
equilibrium electron distribution consists of a Maxwellian part with 
fcTc = 93 keV and a power-law tail with the slope determined by 
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Figure 1. Equilibrium spectra (upper panel) and electron distributions 
(lower panel) for SSC models with the varying magnetization parameter: 
r(B = 0.1 (black solid line), 1.0 (red dotted line), 10 (blue dashed line). 
Other parameters are fiducial. Runs 1-3 from Table[T] 



the radiative cooling s — Finj + 1 = 4. The synchrotron emission 
coming from thermal electrons is completely self-absorbed and the 
low-energy part of the spectrum arises from the synchrotron pro- 
duced by a non-thermal electron population only (with Lorentz fac- 
tors 7 > 30). On the other hand, the high-energy Comptonization 
component, cutting off at ~ 100 keV, is mostly produced by the 
Maxwellian part of the electron distribution. Above the cut-off en- 
ergy, there is a high-energy tail resulting from Comptonization on 
the non-thermal electrons. 



3.1.1 Dependence on the magnetic field 

Let us consider the influence of the strength of magnetic field on 
spectral formation. The energy density of typical magnetic fields 
in the vicinity of black holes is assumed to be of the order of 
the equi parti tion with the radiation field or gas energy density 
dOaleev et al.lll979h . In Fig. [T] the equilibrium spectra and elec- 
tron distributions are shown for rye = 0.1, 1.0, 10. When the role 
of the magnetic field increases, it causes the increasing emission in 
the low-energy part of the spectrum, where the main contribution 
comes from the synchrotron mechanism. Subsequently, the electron 
cooling occurs faster (for the comparison with the re-acceleration 
time see equation |9}. This implies the decreasing role of inverse 
Compton scattering on non-thermal electrons and thus leads to a 
weaker high-energy tail in the photon distribution. The resulting 
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Table 1. Model parameters and results. 
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0.1 





0.800 


93 


28 


m 


10^ 


0.03 


10 


3.0 


1.0 


0.1 


0.1 


0.898 


63 


29 


m 


10^ 


0.1 


10 


3.0 


1.0 


0.1 


0.3 


1.030 


36 


30 


m 


10^ 


0.3 


10 


3.0 


1.0 


0.1 


1.0 


1.191 


16 


31 


m 


10^ 


1 


10 


3.0 


1.0 


0.1 


3.0 


1.335 


4 


32 




108 


1.0 


10 


3.5 


1.0 


0.25 


10 


1.640 


3 



" The fiducial parameter set. '' The value of the temperature is not well defined as the electron distribution is 
poorly fitted by the Maxwellian. The spectrum was calculated including bremsstrahlung. 



Comptonized spectrum becomes softer. As electrons lose more en- 
ergy in synchrotron emission and cool more efficiently, the non- 
thermal tail drops and the thermal part of distribution extends to 
higher energies. At the same time, the peak of the thermal distribu- 
tion is determined by Coulomb processes; thus the equilibrium tem- 
perature is independent of the magnetization. Although the electron 
temperatures are almost identical in these simulations, the slopes 
are different. This is due to the fact that the slope is determined by 
both thermal and non-thermal electron populations, with a harder 
spectral index in the case of more powerful electron tail. 

One can notice in Fig.lHb), that the normalization of the elec- 



tron tail at 7 < 10 does not scale inversely proportional to magne- 
tization. The competing processes determining the normalization 
(and shape) of the high-energy tail are Compton and synchrotron 
cooling. Their relative role in tail formation is determined by the 
magnetization parameter (see equation |6]l. With the increasing 
role of the magnetic field, the cooling becomes dominated by syn- 
chrotron processes, and from this moment the normalization scales 
inversely with r\^. Also at high Lorenz factors, Compton scatter- 
ing plays a less important role, because cooling occurs partially in 
Klein-Nishina regime. 
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Figure 2. Equilibrium spectra (differential compactness, upper panel) and 
electron distributions (lower panel) for SSC models with varying masses 
of the black hole: M = 10 (black solid line), 10^ (red dotted line), 10^ 
(blue dashed line), lO'^M© ( green dot-dashed Hne) and 10^ Mq (yellow 
three-dot-dashed line). Other parameters are fiducial. The broadest spec- 
trum corresponds to the highest mass. The magenta line corresponds to the 
power-law of the energy index a = 5/6. Runs 1, 4-7 from Table[T] 



3.1.2 Dependence on the mass 

Variations in the mass correspond to variations of two parameters 
simultaneously: the Schwarzschild radius and the Eddington lumi- 
nosity. If we fix the Eddington ratio and the size (in Schwarzschild 
radii), then the compactness parameter becomes independent of the 
mass. The spectral energy distribution Ie is not expected to change 
much. However, the temperature of the accretion disc photons as 
well as the magnetic field strength (and characteristic synchrotron 
frequencies) do depends on mass. 

If the seed photons for Comptonization come from an accre- 
tion disc, the size and disc temperature scale as 7? oc Mbh and 
Td oc A/gjJ''*, respectively (see, e.g. iFrank, King & Rainell2002h . 
For the black holes with higher masses and the same compactness, 
the Comptonized spectrum becomes broader and the X-ray spectra 
of Seyf erts are expected to be softer than those of GBHs at the same 
ratio / jZdziarski et"ai]|2003h for Q < 1. 

The situation is more complicated for the SSC mechanism. 
The low-energy part of the spectrum is determined by synchrotron 
radiation which is strongly self-absorbed up to some critical fre- 
quency, the so-called turnover frequency i/t- For hybrid electron 
distribution, containing Maxwellian and power-law parts, there are 
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Figure 3. Same as in Fig.|2ja), but for tjb = 1-0. Runs 8-12 from Table[T] 



lMahadevarj|l997l: IWardzinski & Zdziarski l l2000l) and power- law 
(e.g. lRvbicki & Li ghtmanll 19791 : Wardziiis ki & Zdziarskill200lh 



lyf oc B^tI oc oc Mj^^'+^ , (13) 

where we used equation l[2]l, obtaining the last scaling. Here s = 
Tinj -I- 1, if the cooling of particles emitting around the turnover 
frequency is dominated by synchrotron radiation or Compton scat- 
tering and s — Pinj — 1, if Coulomb collisions dominate (see equa- 
tionll2t. In our simulations, the synchrotron radiation coming from 
the Maxwellian part of the electron distribution is completely self- 
absorbed and only non-thermal electrons are responsible for the 
synchrotron radiation coming from the region. Using equation lO 
we then deduce (for Pinj =3.0) 

0.33H-0.38 



l^t OC 



R 



(14) 



two corresponding turnover frequencies: thermal ul^ oc b'^'^^ (e.g. 



The lower and upper limits of the power-law index correspond 
to Coulomb and radiative regimes, respectively. Moreover, for the 
pure SSC simulations, the synchrotron-emitting particles are in the 
radiative cooling regime (i.e. the cooling by Coulomb is negligi- 
ble). The turnover frequency is generally smaller for more massive 
objects (see Figs|2}0 and lower magnetizations. The latter effect 
can be seen in Fig.[T] however, the dependence on rjB is weak since 
the effects of the increasing magnetic field and decreasing normal- 
ization of the tail of the electron distribution due to stronger cool- 
ing nearly cancel each other out. We find that for different values 
of magnetization parameter, the slope of the thermal Comptoniza- 
tion component can become harder or softer with increasing mass 
(Figs. |2] and [3]l, depending on its initial value and the slope of the 
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Figure 4. Same as in Fig.O but for L/L-^dd = 10~^ and r = 0.3. Runs 
13-17 from Table[T] 

non-thermal electron distribution. The electron distribution is inde- 
pendent of the mass due to the fact that cooling and heating (re- 
acceleration) times scale with the mass similarly, thus their ratio is 
independent of the mass (see equations l9land[7ot. 

The synchrotron luminosity at vt can be estimated by remem- 
bering that the intensity at the turnover frequency is just a source 
function of the power-law electrons and that the relativistic elec- 
tron of the Lorentz factor 7 emits most of the synchrotron power at 
frequency u ~ 27^z/b (where ub ~ eB /2-KmcC is the cyclotron 
frequency). Then, th e compactness of the synchr otron photon spec- 
trum at vt scales as tevbicki & Lightm^]! 19791) 

Rv'^'^ s-3 s-3 s-3 

, synch OC OC 

Eliminating the mass using equation ( 113b . we get the dependence 
of differential compactness (at a turnover frequency) on the photon 
energy 

£^^E,aynch OC ^^^tI. (16) 

This corresponds to a power-law of the energy index a = 5/(s+2). 
In Figs|2]and[3l the magenta line corresponds to s = Finj + 1 = 4, 
giving -EZe, synch OC E^^^ (i.e. a = 5/6). For an initial spectral 
index a, increasing the mass of the compact object leads to softer 
SSC spectra if s < 5/a ~ 2. For the canonical value a « 0.7 
in the hard-state spectra of GBHs we get s < 5 (or equivalently 
Tinj < 4). For a steep spectral slope to start with at 10 M© (for 
instance, the slope as in Fig.|3), we get a more stringent constraint 
then s < 5/0.908 - 2 3.5 and Vinj < 2.5. Thus, a hard elec- 




log E (kcV) 
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Figure 5. Equilibrium spectra (upper panel) and electron distributions 
(lower panel) for SSC models with varying optical depths r = 0. 1, 0.3, 1, 3, 
10. The hardest spectrum and the coldest electron distribution correspond 
to the highest t. Other parameters are fiducial. Runs 1, 18-21 from Table[T] 



tron injection function is required to get the systematically softer 
spectra in Seyferts in terms of the SSC mechanism, keeping all the 
parameters (except Mbh) constant. 

We emphasize here that whether the spectrum hardens or 
softens with the increasing mass depends both on the initial (for 
A/bh ~ IOM0) slope and the scaling of the seed photon compact- 
ness. The power-law slope of the Comptonization component will 
tend asymptotically to the value a — 5/(s + 2) (thus, hardening 
or steepening) with the increasing mass of the object independently 
of other parameters, leading to a stability of the spectrum. In mod- 
els with the seed photons coming from the accretion disc (with / = 
const) this asymptotic slope is a = 1. 

In Fig. |4] we plot spectra and electron distributions obtained 
for the same parameters as in Fig. [2] but for lower r and L/Lsdd. 
The tail in the electron distribution becomes relatively stronger re- 
sulting in more luminous synchrotron emission. This softens the 
Comptonized part of the spectrum (compared to the cases in Fig.|2l(. 
The slopes of electron distributions between momenta p — I and 
10 resemble the injection function, implying inefficient cooling at 
these energies. Below the peak, the Coulomb collisions dominate 
the re-acceleration, leading to a formation of the Maxwellian dis- 
tribution, and for p > 10, radiative cooling softens the electron 
distribution. 
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Figure 6. Dependence of the spectral slope on the optical depth for various 
sets of parameters (Mbh/Mq, L/LEdd. rinj)= (10'', 10~^, 3.0) (solid 
line), (10^, 10"^, 3.0) (dashed line), (lO'', lO'^, 2.0) (dot-dashed line). 
The corresponding dependence for GBHs with parameters (10, 10~^, 3.0) 
is shown with the dotted hne. Other parameters are fiducial. 
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Figure 7. Dependence of the spectral slope on the Eddington ratio for for 
pure SSC models, (a) SMBH (Mbh = 1O^M0, r?B = 1.0) and (b) GBH 
(^BH = IOMq, ?7b = 1.0). Different dependencies of optical depth on 
the luminosity are shown: t oc L^/^ (soHd line), r oc L (dashed line), 
T =constant (dotted line). The curves intersect at L = 10~^LEdd where 
T = 1.0. The light (green) curves correspond to all radiation processes 
taken into account, while for the black curves the bremsstrahlung is ne- 
glected. 
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Figure 8. Equilibrium spectra (normalized to the total compactness, upper 
panel) and electron distributions (lower panel) for pure SSC models with 
varying Eddington ratio: L/Z/Edd = 10~^ (black solid line), 10~^ (red 
dotted line), 10^^ (blue dashed line). The optical depth is constant in these 
simulations. rjB = 1-0, other parameters are fiducial. Runs 2, 22 and 23 
from Table[T] 



3.1.3 Dependence on the optical depth 

Now let us consider the role of the optical depth. Again, we take 
the fiducial parameters and let r vary. As the optical depth grows, 
the electrons have more time to cool before being re-accelerated, 
leading to a drop in the electron temperature (Fig.[5}3). Lower tem- 
perature and higher optical depth lead to efficient thermalization 
by Coulomb collisions, driving a larger fraction of particles to a 
Maxwellian distribution. As the electron temperature decreases, the 
peak of the thermally Comptonized spectrum also moves to lower 
energies (see Fig.|5^). At the same time, the number of particles in 
the non-thermal tail of the electron distribution decreases, leading 
to weaker partially self-absorbed synchrotron emission. The source 
therefore becomes photon- starved and the spectrum hardens. 

At low optical depths the particles do not have time to cool 
between re-accelerations and the electron distribution resembles 
the injection function with a hard tail. The synchrotron emission 
is therefore strong and dominates the emission in the 0.1-10 eV 
range for r = 0.1. 

The effect of the hardening of the spectrum with increasing 
optical depth was found for both GBHs and SMBHs irrespective 
of the values of the other parameters (see Fig.[6]l, suggesting that it 
has the strongest influence on the slope. 
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Figure 9. Equilibrium spectra (upper panel) and electron distributions 
(lower panel) for pure SSC models with varying injection slopes: Fi^j = 
2.0 (red dotted line), 3.0 (black solid line), 3.5 (blue dashed line). Other 
parameters are fiducial. Runs 1, 24 and 25 from Table[T] 



3.1.4 Dependence on the Eddington ratio 

The change in the Eddington ratio is closely connected to the 
change of the accretion rate, which, in turn, affects the optical 
depth. The latter is rather obvious: if the accretion time is con- 
stant, the higher is the matter supply, the higher is the density and 
therefore the optical depth. We assume that r cx Ad. For different 
types of accretion modes (corresponding to efficient and inefficient 
flows), the observed X-ray luminosity depends on the accretion rate 
in a different manner. To investigate spectral formation at various 
accretion rates, we consider a simple power-law scaling 



To 



L 



(17) 



where ro = 1.0 and 6 = 0.5, which corresponds to radiatively in- 
efficient accretion (r oc M and L oc AI^), 9 — 1.0 (corresponding 
to radiatively efficient accretion) and ^ = 0, i.e. constant optical 
depth (ad hoc). 

The results of simulations are shown in Figs|7]and[8] The spec- 
tral index is almost independent of the Eddington ratio, if the op- 
tical depth remains constant and bremsstrahlung is not accounted 
for The turnover frequency is smaller for lower Eddington ratios: 
the fraction of the magnetic energy (77B) remains constant; there- 
fore, the magnetic field itself decreases as i? cx (L/LEdd)^'^^, 
leading to the reduction of vt. The spectral slopes in the 2-10 keV 



range are determined by Comptonization on thermal electrons. Our 
simulations show that for the case of constant optical depth, the 
Maxwellian part of the distribution remains almost unchanged. The 
Comptonized slope can be approximated as (Beloborodov. 1999a): 



9 _ 



2/9 



(18) 



where y — 4(0c + 40^)r(T + 1) is the Compton parameter and 
Oo — kTc/mc^ is the dimensionless temperature of electron dis- 
tribution. Hence, for constant optical depth, equal temperatures of 
the Maxwellian distribution give the same slope of the Comptoniza- 
tion spectrum (see Fig. (8)1. When the optical depth changes together 
with the Eddington ratio, the electron temperature is no longer con- 
served, and the spectral slope changes too. At sufficiently low lu- 
minosities and relatively high optical depths the X-ray spectrum is 
no longer dominated by Comptonization, but instead determined 
by bremsstrahlung. 

Because the injection and cooling rates scale the same way 
with the Eddington ratio, the electron distribution at energies pro- 
ducing synchrotron emission above the turnover frequency 



4.2 X 10" 



(19) 



has a power-law shape with s = Finj + 1 and does not change in the 
cases plotted in Fig. [8] Using equations il3l and J15t , one gets the 
dependence of EIb/I (the spectral distribution scaled to the total 
compactness) on the turnover frequency for the constant R: 



EIe RpV"^ 



(20) 



giving a = 5/(s + 2) = 5/6, which happens to be very close to 
the X-ray spectral index 0.89. Thus, the spectral slope in this case 
is conserved. 

The three electron distributions differ much in the range of 
Lorentz factors between 7 ~ 3 and 20. The effect is due to the in- 
creasing role of Coulomb cooling for lower Eddington ratios (com- 
paring to radiative cooling, see equation[8}. For the adopted param- 
eters and 7 « 10, the ratio of radiative (Compton) and Coulomb 
time-scales is 

^cooljCoul 



tc 



2.8Zinj. (21) 



Thus, for L/Lsdd ~ 10 
inates, but for L/L^dd = 



^ (lini ~ 10), the Compton cooling dom- 
10"* (with hnj ~ 0.1), the slope of the 
electron tail at 7 between 3 and 20 is predominantly determined 
by Coulomb collisions and s « Finj — 1 = 2 (see equation [T2J. 
Coulomb collisions also cause the heating of thermal electrons to 
slightly higher temperatures. 



3.1.5 Dependence on the injection slope 

The exact mechanism responsible for the particle acceleration in 
the inner parts of accretion flow is poorly understood. It can be 
related to weak shocks within the accretion flow or reconnection 
of the magnetic field. The slope of the injected electrons is strongly 
model-dependent. On the observat ional side, the slop e in GBHs can 
be estimated from the MeV tails jMcConne ll et al."2002'). but for 
AGNs the observations are not yet very constraining ( Gonde k^t al.l 
1 1994 1 Johnson et aljl99'1) . 

To investigate the dependence of the results on the injection 
slope, we consider three cases with Finj = 2.0, 3.0 and 3.5. The 
variable injection slope influences the spectrum mainly through its 
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Figure 10. Spectra of a SMBH, where the size of the active region is as- 
sumed to be _R = 300-Rs ■ Runs 26 and 27 from Table [T] Spectra ai'e cal- 
culated without accounting for bremsstrahlung processes (black solid line) 
and including those (blue dotted line). 

strong effect on the strength of partially self-absorbed synchrotron 
radiation produced by non-thermal electrons. Generally, the syn- 
chrotron spectrum fro m the hybrid electron distribut ion can be di- 
vided into four parts dWardzinski & ZdziarskilllOOlh : at low ener- 
gies, emission and absorption are dominated by thermal electrons 
(Rayleigh-Jeans spectrum), at slightly higher energies, the emis- 
sion from non-thermal electrons becomes dominant, but the ab- 
sorption is still determined by the thermal population. Then the 
non-thermal absorption becomes more important, and the spectrum 
is proportional to the source function of non-thermal electrons. Fi- 
nally, at high energies, the non-thermal emission becomes optically 
thin with the spectral index a = [a — l)/2. 

The cases Finj — 2.0 and 3.5 are qualitatively different, as in 
the former case the number of electrons emitting above the turnover 
frequency is large and the synchrotron spectrum is fiat (with a = 1 
corresponding to s = 3, see Fig.|9}. This leads to a large amount 
of soft photons available for Comptonization, strong cooling of the 
electrons and a softer Comptonized spectrum. For high values of 
Finj , the number of electrons emitting above the turnover frequency 
is reduced together with the synchrotron luminosity causing higher 
electron temperature and a harder Comptonized spectrum. 

3.2 The role of bremsstrahlung in SSC models 

We have discussed the spectral formation under the joint action 
of synchrotron and Compton processes. Let us now consider how 
the free-free emission affects the resulting spectra. Major contribu- 
tion comes from electron-electron and electron-proton collisions 
(the amount of positrons is negligibly small in our simulations; 
thus, the number of protons is approximately equal to number of 
electrons). The bremsstrahlung emission coefficients depend on the 
square of the number density of electrons, while the Compton emis- 
sivity depends on the product of the number density of particles 
and photons. Generally one would expect the former processes to 
play dominant role in the medium with high particle and low pho- 
ton number densities. The relevant scalings are oc t/R, and 
riph oc Ltcsc/R^ oc L/B? . Then, the bremsstrahlung emission 
coefficient scales as jbr oc /R^, and Compton scattering emis- 
sion as jcs oc Lt/B^ . The relative importance of bremsstrahlung 
emission is thus characterized by the ratio 
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Figure 11. Equilibrium spectra and electron distributions for different Ed- 
dington ratios (from the top to bottom) L/LEdd = 0.01, 0.03, 0.1, 0.3 and 
1 , and varying disc contribution / = 0, 0.1,0.3, landS.Runs l,28-31from 
Table [T] The disc temperature is assumed to be = lO** K. 



jbr/jcs OC tR/L oc t/L (22) 

The effect of bremsstrahlung on the spectral slope for different de- 
pendencies of T on -L/I/Edd can be seen in Fig. |7] From equa- 
tion l l22t we see that for a given Eddington ratio (and size), its im- 
portance primarily depends on the optical thickness of the source. 
In the cases considered in Fig.]?] the free-free radiation only affects 
the slope in the case when r is kept constant (i.e. large). Thus, 
bremsstrahlung processes can be energetically important for the 
system of larger size (and small compactness). This can correspond 
to some type of radiatively inefficient accretion jNaravan et al.l 
[T998). An example of the resulting spectrum is shown in Fig. 1101 
If the optical depth decreases together with the Eddington ratio, 
bremsstrahlung emission remains negligible. 

However, even if bremsstrahlung is energetically not signifi- 
cant, it still can produce a hardening of the otherwise power-law- 
like spectrum above a few keV making the overall spectrum con- 
cave. Such a spectral complexity might ha ve been a reason for ex - 
istence of soft excess in a GBH Cyg X-1 llbragimov et alj|2005[) . 
In Seyferts, the data are not of the same quality to detect such a 
hardening. 
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3.3 Influence of the accretion disc and slope - luminosity 
correlation 



OBSERVED SPECTRAL PROPERTIES OF 
SUPERMASSIVE BLACK HOLES 



So far have considered pure SSC models, where the only source 
of soft seed photons is synchrotron radiation. This regime can 
be realized at sufficiently low accretion rates, when the disc is 
far away from the region in which the majority of energy is lib- 
erated. For black holes in X-ray binaries, there is a strong evi- 
dence that the disc inner radius changes with the accretion rate. 
The spectral slope correlates with the amplitude of reflection, the 
wid th of the iron line, and the characteristic frequencies of variabil- 
ity jZdziarski et alJll99^;lGilfanov Churazov & Revnivtsevll 19991 : 



iRevnivtsev. Gilfanov & Churazovl 120011 : llbragimov et alj l2005h . 
Similar correlations are found in Seyfert ga laxies in individual 



sources jMagdziMZ et alj '1998: Nandra et al. 2000) as well as m 
the set of^GNs jZdziarski et aL1999.:.Lubinski & Zdziarskil200ll : 



lMatlll200ll : rPapadakis et al.ll2002l : lzdziarski et al.ll2003l) . 

The spectral transitions in GBHs are characterized by the dra- 
matic change in the spectral shape with corresponding changes 
in the electron distribution. In the hard state, the electrons are 
mostly thermal probably only with a we ak tail which is respon - 
sible for the production of MeV photons jMcConnell et al.ll2002h . 
In the soft state, the electrons have mostly non-thermal, power- 
law-like distribution prod u cing long tail ex tending up to 10 MeV 
jPoutanen & Coppil 1 19981 : iPoutane j Il998h . The strong variation 
in the electron distribution can be caused mostly by the increas- 
ing Compton co oling when the disc moves in as the accretion 
rate increa ses ( Vurm & Poutaiienl l2008l : iPoutanen & VumJ 1 20091 : 
iMalzac &"Belmon t 2009). Here we repeat simulations presented in 
aforementioned papers, but for the parameters typical for SMBHs. 

Once the disc moves sufficiently close to the X-ray emitting 
region, the seed photons from the disc start dominating over the 
internally produced synchrotron photons. We simulate this pro- 
cess by injecting additional blackbody photons to the simulation 
volume. We parametrize the importance of these photons by the 
relative fraction to the luminosity dissipated in the active region 
/ = Ldisc/iiiij. These seed photons are injected homogeneously 
into the region. We expect that / is a strong function of the cold 
inner disc radius, which in turn depends on the accretion rate and 
total luminosity. We assume a power-law dependence 



(23) 



where /? and fo are the parameters to be determined from the data 
or theory. The results of simulations for /3 = 1 and fo — 3 are 
shown in Fig.[TT] As the fraction / increases, the disc quickly re- 
places synchrotron radiation as a dominant source of seed photons 
for Compton upscattering. Strong soft radiation field implies faster 
cooling, leading to a decrease of the electron temperature and a 
more pronounced high-energy tail. This, in turn, leads to the de- 
cline of the Compton parameter associated with the thermal popu- 
lation and the spectrum softens (see equation (Ts}. The role of ther- 
mal Comptonization in the formation of the high-energy spectrum 
becomes progressively weaker and is taken over by Compton scat- 
tering on the non-thermal population of electrons. For high /, the 
spectrum is dominated by the disc blackbody, whereas the soft non- 
thermal tail extends to several MeV without a cut-off. The shape of 
the spectrum above ~ 100 keV strongly depends on the assumed 
maximum electron Lorentz factor 7ho. 



Among first surveys on X-ray spectral indices of AGNs, 
i Nandra & Pounds! ( 1 1994!) have noticed that majority of Seyfert 1 
galaxies have spectral energy index a = 0.95 with small disper- 
sion, a = 0.15. Ho wever, addition of NLSyls bro adened the range 
of spectral indices jBrandt, Mathur & Elvi3ll99"^ . After that, the 
energy indices a ^ 1.0 were also f o und for z < 6 type I ra dio-quiet 
quasars (e.g.. IPorguet et al.l l2004l : IShemmer et alj hoO^. Lately, 
it became clear that dispersion depends on the sample: the more 
o bjects are considered , the wider is the spread (see, e.g.. Fig. 8 
in iBurlon et l|). Also recent studies revealed tight corre- 

lation between the spectral sl ope and the luminosity of the ob 
ject in Eddington units (e.g., Wang, Watarai & Mineshige 2004|; 



Shemmer et"aLll2006l: FVasudevan & Fabjan, 2001; Middlgton e^al 



2008!: ISobolewska & PapadakisI 12009! : IVasudevan & Fabian|[2009l: 
Zhou&ZhadbOlO!) . In contrast, the low-luminosity AGNs (e.g. 



LINERs), tend to have rather stable spectral slopes with the aver- 
age Q ~ 1. 

A standard view is that at low mass accretion rates the vicin- 
ity of the bl ack hole is occupied by a hot radiatively inefficient ac- 
cretion flow jNaravan et al.ll998l : miarr Ouataert & Naravarj200"3l: 
[Yuan 2010). In the case of very large accretion rates typical for 
NLSls and quasars, the standard cold accretion disc extends all 
the way down to the last stable orbit around the central black hole 
with the hard X-rays being probably produced in the accretion disc 
corona. In the intermediate case represented by Seyfert galaxies, 
the cold disc is probably truncated and the X-rays are produced in 
the inner hot flow. Let us now apply our model to the observed 
properties of AGNs at different accretion rates. 

We assume the dependence f{L) given by equation l |23t to- 
gether with t{L) given by equation l |17| l. The coefficients and 
the indices can be obtained from the direct comparison with the 
data. Fig. [12] (solid line) shows the dependence of the spectral 
slope on the Eddington ratio for a model with r — {L/L^dd)^^'^ 
and / = 2(L/LEdd). The model describes well the observed 
slopes for a very broad range of luminosities, as we show be- 
low. The self-absorption frequency is sufficiently low for objects 
with high masses and the synchrotron emission falls in the infrared 
range. This model predicts an intimate relation between the in- 
frared/optical and X-ray variability in low-luminosity objects. 



4.1 LINERs and stability of the spectral slope 

Low-ionization nuclear emission-line regions (LINER s) are de- 
fined by narrow optical emission lines of low ionization jHeckmanl 
|l980). Some fraction of LINERs are believed to host an AGN, 
with mass accretion rates far below Eddington. The X-ray spec- 
trum can well be fitted with a s ingle power-law of a ~ 1.0 — 1.2 
jGonzalez-Martm et al. I l2009d) , which shows very little intrin- 
sic absorption and weak (or absent) signatures of Fe Ka emis- 
sion or Compton reflection (Ho 2008, and references therein). 
Thus, the data suggest no accretion disc in the vicinity of the 
central black hole, making LINERs appropriate candidates for 

pure SSC mechanism. We have taken the spectral indices from 

I ' ^ ' 

the sample of IGonzalez-M artm et al. ( 2009b) and the correspond- 
ing Eddington ratios from Satyapal et al. ( 2005) (where avail- 
able). Additimalyj_^o£_sonie objects, we took the black hole mass 
from [ Gonzalez-Martm et al] ll2009£l) and X-ray luminosity from 
iGonza lez-Martm et al. ( 2009b!) and calculated Eddington ratios as- 
suming typical LINER X-ray-to-bolometric correction k = 1.8. 




Figure 12. Dependence o f the spectral slope on the Eddington ratio. The red diamonds mark LINERs (compi lation from ISatvapal et alj |2005| ; 
iGonzalez-Martm et al.| |20o"9 b). Other sy mbols correspond to the dat a on Seyfert galaxies and quasars: magenta crosses jZhou & Zhadl201ol) , blue squares 
iZhou & Zhang 2010) and green triangles iVasudevan & Fabianl2009h . offset up by 0.2 to coiTect for reflection. The black filled circles show the averaged data 
(binned into intervals of one decade in L/L^^^) with the con'csponding dispersion. The black solid line shows the result of our simulations with parameters 
/o = 2 and /3 = 1 from equation (23) and tq = 1.0 and 9 = 0.5 from equation 1171 . The dashed line corresponds to /o = 3 and /3 = 1 and the constant 
optical depth r = 1.0 (corresponding spectra and electron distributions are shown in Fig. II It . 



Then, the data were binned into sections of one decade in L/ L^dd 
and averaged. Individual data points as well as the average spectral 
indices (and their dispersion) are shown in Fig. [T2] The average 
spectral index is consistent with a stable value ~ 1.0 and shows no 
correlation with luminosity. 

Our simulations show that the spectral slopes produced by a 
SSC mechanism are very stable. Among all the parameters, the op- 
tical depth affects the slope most dramatically. For a wide range 
of parameters changing by orders of magnitude {rjB ~ 0.1 — 10, 
r = 0.1 — 3.0, F = 2 — 4), the corresponding spectral index vari- 
ations are within Aa < 0.4. The typical values are ~ 0.8-1.0 and 
do not depend very much on the mass of the central object. We also 
found that increasing the central mass at the same Eddington ratio 
leads to the asymptotic index a = 5/(s + 2). For the steady-state 
index of the electron distribution s = 3 — 5, we have a ~ 0.7 — 1, 
again consistent with the observed hard spectra. This argues in fa- 
vor of the SSC mechanism as the origin of the broad-band spectra 
in LINERs. 



4.2 Seyfert galaxies 

A typical X-ray spectrum of a Seyfert galaxy can be decomposed 
into an underlying power-law and a component arising from re- 
flection and reprocessing of the intrinsic emission. We are inter- 
ested in formation of the intrinsic emission. The low-luminosity 
end of Seyferts distribution at L/L-^dd ~ 10~^ has intrinsic spec- 
tral slopes a ~ 1 very similar to that of LINERs (see Fig.ll2t. Thus 
these sources can also well be described by pure SSC model. 

We note, however, that Seyferts demonstrate spectral vari- 
ations with varying flux. The dependence between the X- 
ray spectral slope and the X-ray flux fo r individual objects 



can often be fitted 
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with a power-law 
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_ Lamer et al.l l2003h . Similar correlations 

were found for GBHs, for example, in Cyg X-I jZdziarski et al.l 



I2OO2I) . If we interpret the observed correlation in terms of the vary- 
ing disc contribution ./(-L) given by equation i23t and assume that 
the spectral slope is determined by the accretion rate only (within 
small uncertainties corresponding to specific parameters of the sys- 
tem), the data require /3 > 2. This implies a strong dependence of 
the inner radius of the accretion disc on the mass accretion rate. 



More recently, ISobolewska & PapadakisI ( |2009|) presented the 
statistical analysis of the correlation for ten AGNs. The average 
a — F correlation was shown to be less steep than that for each 
individual object. This effect might be caused by the fact that the 
same spectral slope in their data correspond to different Edding- 
ton ratios in different objects. Thus tight correlation for any indi- 
vidual object is somewhat smeared out when considering an en- 
semble of sources. The reason of this smearing might be the badly 
determined mass of the central object. A similar slope-luminosity 
(in Eddington units) cor relation is also observed in a large sample 
of Seyferts and quasa rs 1 Vasudevan & Fabian 2009 : Zhou & Zhaol 
I2OICI; IZhou & Zhand l2O10h as shown in Fig. [12] We can inter- 
pret this correlation in terms of the varying contribution from the 
disc to the total luminosity and an increasing role of the disc pho- 
tons for Comptonization. The dependence of the spectral slope on 
the Eddington ratio is shown in Fig. [12] Two cases are consid- 
ered: (1) varying optical depth r — {L/ L^dd)^'''^ with the disc 
luminosity fraction / — 2{L/L-Edd)\ and (2) the constant opti- 
cal depth r = 1.0 with / = 3 (L/I/Edd) (the spectra and elec- 
tron distributions for this case are shown in Fig.ll It. In both cases 
strong correlation a — L/L^dd wa s found. The best-fit ting rela- 
tion a oc 0.3 log (L/I/Edd) found bv lZhou & Zhaol J2010h matches 
well the corresponding lines from our simulations. Thus at high lu- 
minosities L/I/Edd Si 3 X 10^^ the observed spectral softening 
can be explained by the motion of the inner radius of the accretion 
disc towards the black hole and the increasing role of the accre- 
tion disc in supplying soft seed photons for Comptonization in the 
X-ray emitting region (corona or inner hot flow). 
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Figure 13. A hybrid Comptonization model for NLSyls and quasars. The 
sohd hne corresponds to the total spectrum. Contribution from different pro- 
cesses: synchrotron (dotted line), Compton scattering (dot-dashed line, sum 
of thermal and non-thermal), disc black-body (dashed line). The parameters 
are Mbh = 10*Mq, L/L^dd = 1-0, r = 1.0, ??b = 0.25, Tinj = 3.5 
and / = 10. 



4.3 Narrow line Seyfert 1 galaxies and quasars 

Narrow line Seyfert 1 galaxies seem to form a distinct class of 
Seyferts. They are identified by narrow optical permitted and for- 
bidden lines. The UV properties are unusual for presence of both 
low- and high-ionization lines. Th eir X-ray spectra are on aver- 
age s teeper than those of Seyfert Is ( IMiddleton. Done & Gierliiiskil 
120071) . with commonly observed strong soft excess and frequently 
detected high-amplitude variability. The Eddi ngton ratios are o n av- 
erage higher than in broad-line Seyfert Is (e.g. lKomossal20()8l) . The 
power spectral density (PSD) spectrum of the NLSyl NGC 4051 is 
similar to PSD of soft state of Cyg X-1, but shifted to lower fre- 
quencies jMcHardv et"ai]|2004 . Similar to the soft state of GBHs, 
one can expect that the two Seyfert classes differ by the presence 
of the accretion disc close to the compact object. If the disc ex- 
tends very close to the SMBH, most of the energy is radiated in 
the UV range and completely ionizes the medium around. Hence, 
the lines in NLSyls can be produced only in the region far away 
from the central ob ject, where the gas velocities are relatively low 
( lBrandtetal]|l994l) . 

Soft excess is an ubiquitously observed feature in the 0.4- 
1 keV range in NLSyls and radio-quiet qua sars. It can be mod- 
elled by cool Com ptonization component jPounds etai]|l995l : 
IVaughan et al.l2002h : however its position in the spectrum seems to 
be independent on the central black hol e mass, temperature ofac- 
creting disc and luminosity of the obiect. ^ierlinski & Done! l [2004l) 
proposed an alternative formation scenario where the excess arises 
from the smeared absorption of highly ionized oxygen and iron, 
present in the hot wind from the accretion disc. For the quasar 
PG 1211-1-143, they obtained the slope of the intrinsic spectrum of 
a « 1.7 in the range 0.4-10 keV. This spectrum can be well repro- 
duced by our model where the luminosity of the seed photons from 
the disc exceeds by ten times the power injected to electrons (see 
Fig. |13t . This results in a soft, power-law-like spectrum in the X- 
ray band produced by non-thermal Comptonization. In this model, 
no cut-off is expected up to high 100 keV) energies. 



5 CONCLUSIONS 

We have studied the spectral formation in hot non-thermal plasmas 
under the conditions relevant to the vicinities of SMBHs in AGNs, 
using a self-consistent SSC model. We show that the SSC model 
can reproduce the exponential cut-off at energies above ~ 100 keV 
and the power-law X-ray continuum with energy spectral indices 
a ~ 0.8 — 1.0 observed in low-luminosity Seyferts and LINERs. 
The entire infrared/optical to X-ray spectrum in this objects can be 
produced by the SSC mechanism, suggesting a strong correlation 
between the two energy bands. No specific conditions are required 
to stabilize the slope for a wide range of parameters (magnetiza- 
tion, optical depth, injection slope, Eddington ratio and black hole 
mass). 

We have demonstrated the difference in scaling of the SSC 
model and the two-phase disc-corona models. While in the lat- 
ter case one would expect systematically softer spectra from ob- 
jects with higher masses, in the case of SSC, there exists a lim- 
iting slope (when the mass of the central object tends to infinity) 
a — 5/{s + 2), which depends only on the slope of the tail of 
the electron distribution. Comptonization spectra can get systemat- 
ically softer for objects of higher masses, if the electron injection 
function is sufficiently hard. 

We have tested the role of bremsstrahlung emission in SSC 
models. We found that even though the process is not energeti- 
cally significant, in certain cases, it leads to the hardening of the 
power-law above a few keV and makes the overall X-ray spectrum 
concave. 

We have also studied the role of additional soft blackbody 
photons from the accretion disc. With the increasing amount of 
soft photons, the equilibrium Maxwellian temperature drops and 
the X-ray spectrum softens. The resulting spectra look similar to 
the GBHs in their soft state, which in the case of SMBHs are likely 
to be represented by NLSyls and quasars. We have found that it 
is possible to explain the spectral slope - flux correlation, widely 
discussed in the literature, by parametrizing the fraction of the disc 
photons in the medium as a power-law function of the Eddington 
ratio. This implies, that the inner radius of the truncated accretion 
disc is a strong function of the accretion rate. 
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